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The possibility that the Fermi scale is the only fundamental energy scale of Nature 
is under serious consideration at present, yet cosmic rays may already have provided 
direct evidence of new physics at a much higher scale. The recent detection of very 
high energy particles with no plausible astrophysical sources suggests that these 
originate from the slow decays of massive particles clustered in the halo of our 
Galaxy. Such particles had in fact been predicted to exist beforehand with mass 
and lifetime in the range required to explain the observations. I discuss recent 
work focussing on experimental tests of this speculative but exciting idea. 

1 Introduction 

The only massive particles in the Standard Model to have survived from the 
Big Bang are nucleons — protons and (bound) neutrons — along with a com- 
mensurate number of electrons to yield the observed charge neutrality of the 
universe]^ Considerations of primordial nucleosynthesis restrict the nucleonic 
contribution to the density parameter to £1n 0.1 and it is widely accepted 
that the dark matter in galaxies and clusters which contributes JIdm ^ 0.3 is 
non-nucleonic and probably composed of a new stable relic particle. There are 
many candidates for the identity of this particle but the most popular notion is 
that it is associated with the new physics beyond the Standard Model necessary 

— 1/2 

to stabilize the hierarchy between the Fermi scale, G F ~ 300 GeV, and the 

Planck scale, G^ 1 ^ 2 ~ 10 19 GeV. In particular theories of (softly broken) low 
energy supersymmetry (SUSY) typically imply that the lightest SUSY partner 
is a neutralino with mass of order the Fermi scale, which is absolutely stable 
if the discrete symmetry termed i?-parity is exactly conserved. Interestingly 
enough the relic abundance of such a weakly interacting particle which was in 
thermal equilibrium in the early universe can account for the dark matter. 

In supergravity theories, there is a new energy scale of O(10 n ) GeV — the 
geometric mean of the Fermi and Planck scales. This is the scale of the 'hid- 
den sector' in which SUSY is broken through gaugino condensation induced 
by a new strong interaction, and communicated to the visible sector through 
gravitational interactions. Following the emergence of superstrings (for which 

a We know now that massive relic neutrinos contribute at least as much as the luminous com- 
ponent of nucleons to the present energy density. However they are unlikely to be the dom- 
inant component of the dark matter, based on arguments concerning structure formation. 
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N = 1 supergravity is the effective field theory) it was realised!!! that the hidden 
sector can .also serve to confine fractionally charged states which are a generic 
prediction □ of string theory. This avoids a serious conflict with the unsuc- 
cessful experimental searches for fractional charges but necessarily implies the 
existence of (integrally charged) bound states with mass of 0(lO n l GeV. In a 
specific construction with ST/ (5) ® U(l) unification, it was notedH that most 
such states would be short-lived but that the lightest such state would only 
decay through non-renormalizable operators of dimension > 8 and thus have 
a lifetime exceeding the age of the universe. This introduces a new candidate 
for the constituent of the dark matter — named "cryptons" — interestingly 
similar to nucleons which too are bound states of fractional charges and can 
only decay through non-renormalizable operators. 

However, just as with nucleons, their cosmological origin is a puzzle. If such 
particles were ever in thermal equilibrium their relic abundance would have 
been excessive since their self-annihilations are rather inefficient. For nucleons 
the problem is just the opposite and their very existence today requires an out- 
of-equilibrium origin. If the same were true of cryptons, their relic abundance 
may well have a cosmologically interesting value .[] It is then interesting to ask 
what the observational signatures of such partinles might be. 

Reviving an old suggestionO, we recognisedu that the most sensitive probe 
would be in extremely high energy cosmic rays (EHECR), specifically in the 
flux of high energy neutrinos which would necessarily be created by crypton 
decays. The best constraint we obtained followed from the upper limit on 
deeply penetrating air showers set by the Fly's Eye atmospheric fluorescence 
experiment; this implied that such particles must have a lifetime exceeding 
~ 10 18 yr if they are an important constituent of the dark matter. As this was 
close to the theoretically expected lifetime in theJ'flipped" SU (5) model, I was 
optimistic enough to suggest in a conference talkQ that "... some improvement 
of these experimental sensitivities can rule out (or detect!) ^uch particles". 

Just a few months later the Fly's Eye array detected U an event, consis- 
tent with a proton primary, but with an energy of (3.0 ± 0.9) XplO 11 GeV. 
This was well above the Greisen-Zatsepin-Kuzmin (GZK) cutoff □ energy of 
~ 5 x 10 10 GeV, beyond which resonant photopion production losses on the 
cosmic microwave background should limit the propagation distance of any 
such strongly interacting particle to less than about a hundred Mpc. Over a 
dozen such events have been detected subsequently by the Akeno airshower 
array (AGASA) as well as HiRes, the successor to Fly's Eye, so the absence of 

b It has recently been noted @ that particles with mass of 0(H- m f) ~ 10 13 GeV — also 
dubbed "wimpzillas" — can be created with a cosmologically interesting abundance through 
quantum vacuum fluctuations during inflation or during the subsequent (re)heating process. 
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the GZK cutoff £3 is now well established. However contrary to the expecta- 
tion that such high energy particles, being essentially undeflected by the weak 
intergalactic magnetic -fields, should point back to their sources, the observed 
distribution on the skyEll is consistent with isotropy. This is quite baffling given 
that that only a few astrophysical sites (active galactic nuclei or the extended 
lobes of radio galaxies) am. capable of accelerating such particles, even in prin- 
ciple, and there are noneEj along the arrival-directions within the propagation 
range. Hence it is generally acknowledged E£l that there is no "conventional" 
astrophysical explanation for the observed EHECR. 

2 EHECR from decaying dark matter 

Faced with the above conundrum, some authors have resorted to desperate 
measures, e.g. postulating that the intergalactic magnetic field may be a thou- 
sand times stronger than usually believed, so capable of isotropising particles 
from a nearby active galaxy. However, following from our previous discussion, 
there is a natural explanation Lj for both the observed isotropy and absence of 
the GZK cutoff if the EHECR originate from the decays of metastable cryp- 
tons which are part of the dark matter.^] This is because such particles will 
behave as cold dark matter (CDM) and hence cluster in the halo of the Milky 
Way with a concentration ~ fO 4 times higher than the cosmic average. The 
local flux of EHECR will thus be dominated by decays of cryptons in the halo, 
implying two distinct observational tests of the hypothesis. First, the energy 
spectrum and cosmpositipa (nucleons, gammas, neutrinos) beyond the GZK 
cutoff will be determined liJ essentially,-by the physics of crypton decays. Sec- 
ond, there will be a small anisotropyl!3 in the arrival directions of EHECR 
since we are located ~ 8 kpc away from the centre of the Galaxy and should 
therefore observe more particles arriving from that direction than from the 
anticentre. There may also be measurable correlations between arrival times 
of high energy nucleons, gammas and neutrinos. 

2.1 Particle candidates 

As noted above, the possibility n af metastable relic particles with mass of 
O(10 n ) GeV had been proposed tJtl befoxe. the observations of EHECR beyond 
the GZK cutoff. An updated discussionta of such particles in string/M-theory 
confirms that cryptons are indeed favoured over other possibilities such as 

c This was independently proposed by Berezinsky, Kachelriefi and iVtlenkin0, without, how- 
ever, a specific particle candidate in mind. Kuzmin and Rubakov|13 also made a qualitative 
suggestion that EHECR may originate from relic particle decays, however they did not make 
the crucial observation that such particles would be highly concentrated in our Galactic halo. 
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the Kaluza-Klein states associated with new compact dimensions (which are 
too short-lived). The most likely candidate is still a neutral pion-like 'tetron' 
composed of four constituents, with a minimum lifetime of 

1 ( M \ 10 

rx^ , 1 

m x \m x J 

where mx ~ 10 12 ~ 13 GeV, and the scale M of suppression of non-renormalizable 
terms is of O(10 18 ) GeV.^| Thus both the mass and lifetime of the candidate 
particle are motinnlBd by topical physical considerations. This is in contrast to 
other proposals EitS where the mass scale is not given any physical motiva- 
tion and the decays are presumed to be mediated by unspecified instanton or 
quantum gravity effects so as to yield a suitably long lifetime. 

2.2 Calculation of decay spectrum 

Nevertheless all such proposals have a common phenomenology in that regard- 
less of the decay mechanism, the spectra of the decay products is essentially 
determined by the physics of QCD fragmentation Ell and has no major astro- 
physical uncertainties. In particular given that the propagation distance in 
the halo is 100 kpc, much shorter than the GZK range of ~ 100 Mpc, the 
EHECR spectrum at Earth will be the same as the decay spectrum (apart from 
the decay photons which will be degraded through scattering on background 
photon fields). Of course the decay mode (e.g. 2-body vs many-body) may 
well play an important role. However in our picture the decaying particle is a 
singlet under Standard Model interactions and has a mass which significantly 
exceeds the Fermi scale so the inclusive spectra of final state nucleons, photons 
and neutrinos should be relatively insensitive to the precise decay channel. 

We can thus imagine that we have say a e + e~ collider at our disposal 
with a centre-of-mass energy i/s sufficient to create a supermassive particle 
such as a crypton, rather than just a Z° as at LEP (Figure [l]). This then 
decays into quarks and gluons which initiate multi-parton cascades through 
gluon bremsstrahlung. These finally hadronize to yield high multiplicity jets 
when the momentum scale of the process drops below Aqcd- In the present 
context we are only interested in the final yields of nucleons, photons and 
neutrinos into which all the produced hadrons will decay. The production of 
different hadron species is quantified by their respective 'total fragmentation 
functions' F h (x,s) — a^da/dx, viz. the probability distributions for their 

d Other ai*thors0 have also considered string candidates for superheavy dark matter, and 
discussed the confinement of fractionally charged particles into baryon-like states in the 
hidden sector and the discrete symmetries required to ensure their longevity. 
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inclusive production as a function of the scaled hadron energy x = 2Eh/\/s. 
These can be factorized as the sum of contributions from different primary 
partons i = u, d, . . . , g: 

—Gtez^WDfofas), (2) 
z 

where Cj are the 'coefficient functions' dependent on the production process, 
and D\ is the 'universal fragmentation function' for parton % — ► hadron j. The 
essential physics in the hard parton cascade is the logarithmic evolution of the 
strorLg_coupling a s with energy and sophisticated techniques have been devel- 
opedEil to handle divergences associated with collinear and soft gluon emission. 
The formation of the final hadrons is however an inherently non-perturbative 
process and can only be described at present bj^ empirical models encoded 
in Monte Carlo event generators, e.g. JETSETE-3 based on the 'string frag- 
mentation' model, or HERW1GE3 based on the 'cluster hadronization' model. 
These also account for the subsequent decay of the hadrons into the observed 
particles, taking into account all experimentally measured branching ratios, 
resonances etc, so can make detailed predictions of measurable quantities. 




Figure 1: The hadronization processed in the decay of a massive particle. 

Although the fragmentation functions are not perturbatively calculable, 
their evolution as a function of the momentum scale is-governed by the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equationE3 

d f 1 dz 

s-q- s D?(x, S )=Y,J —P 3l {z,a s {s))D } t{x/z, S ), (3) 
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where Pji are the 'splitting functions' for the process parton s -> j, Thus 
by measuring the fragmentation functions at one momentum scale, one can 
evaluate them at another scale. As seen in Figure || the DGLAP equations 
predict violations of 'Feynmann scaling' — a softening of the spectrum with 
increasing energy — in good agreement with data, in this case measured at 
PETRA (VS = 22 GeV) and LEP (Vi = 91.2 GeV). 




• i • i < 1 1 1 1 i . 
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Figure 2: Scale dependence of fragmentation functions £j in e+e experiments. 



At small values of x, multiple soft gluon emission gives rise to higher- 
order corrections, which turn out to be resummable by altering the scale in 
the DGLAP equation (0) from s — » z 2 s; this yields a simple Gaussian function 
in the variable £ = ln(l/x): 



xF h (x, s) oc exp 



2a- 



-(£-£ P ) : 



(4) 



which has a characteristic peak at £ p ~ In s/4, with width cr ~ (In s) 3 / 4 . 
Including 'next-to-leading' corrections t£L.these predictions yields the 'modi- 
fied leading log approximation' (MLLA)Ej'EZI which accounts very well for the 
shape of the observed fragmentation functions at small x. In comparing, with 



data one has to further assume 'local parton hadron duality' (LPHD)Ej, viz. 
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that the hadron distribution is simply proportional to the parton distribu- 
tion. Thus the prediction cannot distinguish between the individual hadronic 
species. Moreover although this kinematic region dominates the total multi- 
plicity, it accounts for only a small fraction of the energy in the cascade, hence 
the MLLA spectrum cannot be correctly normalized. 

With this background, we can review what has been done so far to ex= 
plain the EHECR data in terms of decaying halo particles. Berezinsky et ali3 
adopted the gaussian approximation (|J) to MLLA to infer the spectrum of nu- 
cleons from the decay of a particle of mass mx ■ Although this approximation 
is only valid for small x (<C 0.1), these authors nevertheless normalized it by re- 
quiring that J Q dx xF h (x) = /n where /n ~ 0.05 is the assumed fraction of the 
decaying particle mass transferred to nucleons (on the basis of Z° decay data). 
The rest is assumed to go into pions which decay to yield photons and neutri- 
nos, with neutral pions taking a third of the total energy. Figure || shows their 
fit to the EHECR data (multiplied by E 3 for clarity) with decaying particles 
of mass mx — 10 13 GeV which contribute a fraction £x of the CDM density 
in our halo (taken to have a radius of 100 kpc). Their adopted normalization 
then implies a particle lifetime T x /tg = 2 x 10 10 £x, where to ~ 12 x 10 9 yr is 
the age of the universe. Note that the constant suppression with energy of the 
nucleon flux with respect to photons (and neutrinos) follows from the assumed 
proportionality of the nucleon and pion fragmentation functions. 



28 
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Figure 3: Predictedr^hixes from decaying dark matter particles of mass 10 13 GeV according 
to Berezinsky et aZ;|l3 the GZK-suppressed extragalactic proton flux is also shown. 
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Because of the above problems with the MLLA spectrum, we had already 
considered and rejected this convenient approximation and chosen instead to 
embark on a time-consuming calculation of the fragmentation functions in 
the kinematic region of relevance to the data, using the HERWIGE3 event 
generator. In doing so we were initially motivated to test an argument due to 
HillEj that the fragmentation spectrum at large x should be oc (1 — x) 2 . By 
normalizing to the total multiplicity and demanding energy conservation Hill 
was then able to obtain an empirical fragmentation function which could be 
fitted to extant data. Assuming the multiplicity to be oc s 1 / 4 as in the naive 
statistical model of jet fragmentation, this was 

F h = ^x-^{l-xf, (5) 
16 

while by adopting the leading-log QCD prediction for the multiplicity (oc 
exp -\/ln(s/A 2 )), he obtained 

Fh = N{b) eMbVHTMKi-*)\ (6) 

x-\/ln(l/x) 

By fitting this form to PETRA data, Hill found N(b) = 0.08 and b = 2.6. 
On the basis of the same data he assumed that 3% of the hadronic jets form 
nucleons and the other 97% are pions which decay into photons and neutrinos. 




Figure 4: Comparison of the nucleon fragmentation functions fori decaying particles of mass 
10 11 GeV (lowercurve) and 10 3 GeV (upper curve), calculated^ using HERWIG. The Hill 
approximation (tjl) is also shown (dotted line). Note the significant scaling violation. 
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Many authors cJ who have investigated the annihilation of GUT-scalc relic 
topological defects (TD) as the source of EHECR have used these expressions 
to estimate the fluxes. However this is clearly inaccurate since there would be 
large scaling violations (see Figure ||) in going from the PETRA energy scale 
of 22 GeV up_to the very much higher GUT energy scale. This is just what our 
calculationsEJ using HERWIG demonstrate. The functional form (|^) continues 
to provide a good fit to the fragmentation function of nucleons but as shown 
in Figure |J the spectrum becomes significantly softer with increasing particle 
mass, e.g. for mx — 10 13 GjsV, the normalization N(b) drops to 0.0078 (with 
b = 2.8). Thus TD modelso of EHECR which use the fragmentation functions 
(||,^|) overestimate nucleon production by a factor of ~ 10 at high energies. 

The EHECR spectrum in the energy range 10 9 ~ n GeV is well fitted i 
as the sum of two power-laws — the extrapolation of the E~ 3 3 spectrum 
from lower energies and a new flattet, component oc E~ 2 ' 7 which dominates 
above 10 10 GeV. (The AGASA datat3 gwes the slope of the new component 
as — 2.78^tg'33.) There is some indicationcl that the composition also changes 
from iron-group nuclei to protons-at this energy. In Figure [| we see that 
the HERWIG generated spectrum til for a decaying particle mass of 10 12 GeV 
is indeed in reasonable agreement with this new component of cosmic rays. 
Our normalization requires a lifetime Tx/to — 3 x 10 9 £x in the notation of 
Berezinsky et altQ, i.e. a factor of ~ 6 smaller. (This is because they seem to 
have normalized to photons rather than nucleons at 10 10 GeV; see Figure^). 



m v = 10 GeV 
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Figure 5: Expected nucleon fluxQ for decaying halo particle masses of order the hidden 
sector scale compared with data (multiplied by EA for clarity). The spectra are normalized 
at 10 10 GeV to the new component (dashed line)E3 suggested by the observations. 
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However our own calculationstfl suffer from two problems. The first, which 
we were initially unaware of, is that HERWIG has a known tendency Ej to 
overproduce baryons at large x (essentially due to the fragmentation of the 
leading quark from the initial hard process, viz. particle decay in the present 
case) . Although the overall multiplicity is correctly predicted at LEP energies 
(e.g. 0.953 protons per event vs 0.98 ± 0.09 observed in Z° decay), HERWIG 
overproduces nucleons at x *S 0.3 by a factor of ~ 2 — 3. Secondly, in studying 
the evolution of the fragmentation function to very high energies we should 
take into account that the running of the strong coupling a s would be altered 
above ~ 10 3 GeV when SUSY particles begin to be excited from the vacuum. 





Figure 6: Effect of ad hoc baryon suppression factor on HERWIG default (top) 
and isotropic (bottom) baryon cluster decays. 
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Both these issues have been addressed recently in unpublished wor kEl by 
Rubbx. As shown in Figure 6 he finds that agreement of HERWIG with LEP 
datao on baryon production is significantly improved by isotropizing the de- 
cay of the hadron cluster formed from the hard process (rather than have the 
hadron leave the cluster along the direction of the initial quark) . An additional 
suppression of the probability for cluster decay by ~ 20% improves the fit fur- 
ther. To take SUSY particles into account, he evolves the DGLAP equation (||) 
from LEP energies upwards (with the initial sparton fragmentation functions 
calculated with the PYTHIA event generator). The SUSY ft function for 
a s is used, with the flavour thresholds corresponding to the sparticle spectrum 
of a typical minimal supergravity model (with parameters Mo — 800 GeV, 
m 1/2 = 200 GeV, A a = 0, tan/3 = 10, sgn(^i) = +). 
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Figure 7: Comparison of data with spectra predicted using coherent branching: 
Mx = 10 12 GeV without SUSY (top); Al= 10 13 GcV with SUSY (bottom). 



Figure 7 shows his results for both the non-SUSY and SUSY cases. The 
high energy "bump" in our proton spectrum (see Figure has been erased 
but the non-SUSY spectrum continues to reproduce the shape of the data for a 
decaying halo particle mass of O(10 12 ) GeV. The effect of including the effects 
of SUSY on the evolution of the parton cascade is to flatten the spectrum 
further so that a ~ 10 times larger mass is still acceptable. We note that the 
spectral shape differs considerably from the "SUSY-QCD" spectrum calculated 
by Berezinsky and KachelriefiE3 using MLLA. This is not unexpected since as 
emphasized earlier, this approximation is unjustified at large x so cannot be 
normalized (as these authors do) to the energy released in the decay. Moreover 
their assumption of an energy-independent ratio between nucleons and pions 
is invalid; as is evident from Figure 7 this ratio increases with energy. f\ 



3 Conclusions 

Although some progress has been made in sharpening the spectral predictions 
of the decaying halo particle model for EHECR, much work still needs to be 
done. The calculations so far have assumed the simplest decay channel — into 
two partons. However non-renormalizable operators are in fact likely to induce 
many-body decays. The effects of supersymmetry also need to be investigated 
more carefully, e.g. the effects of varying the SUSY parameters and inclusion of 
sparticle decay channels. Nevertheless it is already clear that the general trend 
in the EHECR data can be accounted for by this hypothesis, if the particle 
mass is m x ~ 10 12 ~ 13 GeV and its lifetime is t x ~ 10 16 yr(£x/3 X 1(T 4 ), so 
that even with a very long lifetime such particles need constitute only a tiny 
fraction £x of the halo CDM. It is also clear that TD models E2l, in which mx 
corresponds to the GUT-scale, are already ruled out by the spectral data. 

The next generation of large area cosmic ray, gamma-ray and neutrino 
observatories (Auger, Amanda, Antares, . . . ) is now under construction so it 
is important to refine these calculations in order to make specific predictions 
for the expected fluxes. We emphasize that previous estimates of high energy 
gamma-ray and neutrino fluxes fmm TDc3 are based on the Hill fragmentation 
functions (||]|) pSfhile other worka use the (M)LLA spectrum (§) or its SUSY 
variant. BlasiLJ has calculated in detail the flux of 7-rays in the decaying 
halo particle model but he too uses the Hill and the MLLA spectra. All these 
approximations are inapplicable at the high energies of interest as explained 
earlier, and moreover the spectra of pions are not simply proportional to that 
of nucleons as assumed. Hence it is clear that all these estimates are unreliable. 

e However it never exceeds unity as in our previous work0 using HERWIG which suffered 
from overproduction of hard baryons and gave an incorrect prediction of the p/v ratio. 
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It is essential that further work use the physically more realistic approach to 
calculating fragmentation spectra outlined above in order to devise definitive 
experimental testsu of the decaying particle hypothesis. 
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